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We show that boundaries of 3D weak topological insulators can become gapped by strong interactions while 
preserving all symmetries, leading to Abelian surface topological order. The anomalous nature of the weak 
topological insulators manifests itself in a non-trivial action of symmetries on the quasiparticles; most strikingly, 
translations change the anyon types in a manner impossible in strictly 2D systems with the same symmetry. As 
a further consequence, screw dislocations form non-Abelian defects that trap Z 4 parafermion zero modes. 


INTRODUCTION 


Electronic topological insulators [1-5] display numerous 
exotic properties already at the single-particle level, most fa¬ 
mously protected surface metallicity. Much of the richness in 
these systems emerges from the interplay between symmetry 
and topology. Recently interactions among surface electrons 
have been found to further enlarge the possibilities. In a strong 
topological insulator (STI) the surface spectrum for weakly 
interacting electrons obeying time reversal T and charge con¬ 
servation symmetry features a single Dirac cone. Remark¬ 
ably, strong interactions can fully gap the STI surface without 
violating symmetries [6-9] (as anticipated earlier [10]). The 
symmetrically gapped phases realize non-Abelian topological 
order—the simplest incarnation being the ‘T-Pfaffian’—and 
can be viewed as descending from novel gapless states [11- 
13]. Similar conclusions hold for bosonic topological insu¬ 
lators [14], topological superconductors [15], and topological 
crystalline insulators [16]. (Not all topological systems, how¬ 
ever, admit a symmetric gapped boundary [17].) 

We explore for the first time the fate of strongly correlated 
weak topological insulator (WTI) surfaces. These systems 
comprise ideal settings where one can controllably explore the 
influence of additional symmetries, which enrich the surface 
topological order that we identify in subtle ways and yield an 
interesting interplay with lattice defects. A WTI may con¬ 
veniently be decomposed into a stack of quantum spin Hall 
(QSH) insulators [1-3] with electrons from the helical edges 
tunneling between layers; see Fig. 1. Provided the system pre¬ 
serves time reversal T, charge conservation, and layer transla¬ 
tion symmetry T y , the non-interacting WTI surface hosts two 
massless Dirac cones at distinct momenta [18]. 

Four considerations are useful for anticipating the topolog¬ 
ical order that emerges when interactions gap the WTI sur¬ 
face without violating these symmetries. First, on very gen¬ 
eral grounds the topological order must be anomalous, i.e., 
forbidden in strictly 2D isosymmetric systems. To see this 
consider the thickened torus of WTI depicted in Fig. 2(a), and 
gap the interior surface by interactions but leave the exterior 
gapless. Upon shrinking the torus’s thickness a strictly 2D 
system emerges as in Fig. 2(b). If the gapped surface was 



FIG. 1. (color online) (a) Weak topological insulator built from quan¬ 
tum spin Hall layers. Interlayer tunneling yields two symmetry - 
protected surface Dirac points at momenta (q x ,q y ) = (0,0) an d 
(0, 7 r). (b) Weak topological insulator surface viewed as a bilayer. 


non-anomalous, one could simply strip away the topological 
order, leaving a symmetric 2D system with an ‘impossible’ 
band structure [19] —a contradiction. 

The second consideration regards a domain wall separat¬ 
ing the topologically ordered state from a ferromagnetically 
gapped surface region. The magnetized region carries a non¬ 
zero thermal Hall conductivity and thus the domain wall must 
host gapless modes. In the STI case the thermal Hall con¬ 
ductivity would be half-integer (in units of it 2 k 2 B T /?>h)\ the 
gapless mode’s central charge must also then be half-integer, 
necessitating a non-Abelian topological order. By contrast, 
the two Dirac cones present for the WTI imply an integer cen¬ 
tral charge, suggesting an Abelian minimal topological order. 

The third consideration seemingly contradicts the second 
by raising expectations for a non-Abelian surface state: Since 
interactions can drive the STI’s single Dirac cone into the 
non-Abelian T-Pfaffian [6-9], the WTI’s two surface Dirac 
cones may naturally be gapped independently in the same 
manner. The resulting non-Abelian surface hosts a large num¬ 
ber of quasiparticles—two copies of each T-Pfaffian excita¬ 
tion. As we confirm below, however, the quasiparticle content 
can be reduced via a symmetry-preserving ‘anyon condensa¬ 
tion’ transition that confines the non-Abelian excitations leav¬ 
ing only an Abelian set, consistent with the argument above. 

The fourth consideration results from viewing the WTI as 
a stack of QSH insulators. Any finite stack may be viewed as 
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FIG. 2. (color online) (a) Thickened torus of weak topological insu¬ 
lator with symmetric topologically ordered interior and gapless exte¬ 
rior. (b) Two-dimensional limit where the thickness shrinks to zero. 
The topological order must be anomalous; otherwise one is left with 
an ‘impossible’ 2D band structure. 


two dimensional, with an even-odd effect: the system forms 
a 2D topological insulator with an odd number of layers but 
a trivial 2D insulator otherwise. Since the 2D topological in¬ 
sulator edge cannot be gapped without breaking T or charge 
conservation, this even-odd effect should also appear when in¬ 
teractions gap the stack’s surface to form topological order in 
the limit of infinitely many layers. 

With these points in mind let us investigate gapping of the 
WTI in detail. First we expound on the relation to the STI by 
equivalently describing the WTI surface as a bilayer system, 
partitioning the right- and left-movers from each QSH edge as 
in Fig. 1(b). Here time-reversal T and translations T y inter¬ 
change the two layers while T = TT y does not. Each layer 
represents the surface of an antiferromagnetic topological in¬ 
sulator [20, 21] —which also supports a single Dirac cone— 
and maps to the setup considered in Ref. [11] for studying 
the correlated STI surface. Interactions can drive each layer 
into an electrically insulating, symmetry-preserving ‘compos¬ 
ite Dirac liquid’ that hosts a single Dirac cone built from emer¬ 
gent neutral fermions that carry a fictitious ‘pseudocharge’ e 
(which is not microscopically conserved). Intra-cone Cooper 
pairing of neutral fermions in the bilayer generates ‘double 
T-Pfaffian’ non-Abelian topological order. One can conve¬ 
niently view the resulting anyons as defects in the paired 
condensates; in particular, an he/ 2e vortex carries physical 
charge e/4 and binds a Majorana zero mode. The so-far de¬ 
coupled layers host independent hc/2e vortices. Our anal¬ 
ysis below shows how to obtain an Abelian state for which 
individual h/2 e vortices are confined, yielding e/2 as the 
minimal charge. The Abelian state further includes a neu¬ 
tral quasiparticle—an interlayer ±e/4 dipole—that acquires a 
7 t/ 2 phase when encircling an e/2 excitation. 

We now put this discussion on firmer footing. To facilitate 
gapping the WTI we imagine patterning the surface with 2D 
topologically ordered ‘plates’ as in Ref. [11]. Figure 3(a) de¬ 
picts the decorated structure. Each interface (labeled y) con¬ 
tains a helical QSH mode and two sets of gapless edge states 
from the adjacent plates, one from above and one from below. 
We wish to judiciously select the plates such that (i) local 



FIG. 3. (color online) (a) Weak topological insulator surface dressed 
with 2D topologically ordered ‘plates’, (b) Dragging an any on across 
plates leaves behind ‘invisible’ operators at the interface that get ab¬ 
sorbed into a condensate. The condensates allow the e/2 a quasipar¬ 
ticle to pass freely between plates while the neutral d quasiparticle 
acquires a neutral fermion and thus changes anyon type, (c) Setup 
for discussing weak symmetry breaking. 


interactions within a given interface can remove all gapless 
modes without breaking any symmetries and (ii) the surface 
topological order with minimal degeneracy on a torus appears. 
Note that time-reversal symmetry constrains the latter degen¬ 
eracy to be the square of an integer [22] . 

The interfaces to be gapped are described by a K-matrix K 
and charge vector Q, which specify the statistics and charges 
of low-energy fields, along with a vector X that distinguishes 
Kramers singlets from doublets (for details see [22]). More 
precisely, we have 
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where the ‘ft’ and ‘p’ subscripts indicate quantities for the 
helical QSH modes and plates, respectively. For the QSH 
sector Kh = cr z (here and below a a denote Pauli matri¬ 
ces), qh = (1,1), and Xh = (0,1). For the plates, time 
reversal demands an even-dimensional K p . We assume the 
smallest two-dimensional K p , which can be either fermionic 
or bosonic. We focus on the latter since we find that the 
fermionic case does not permit time-reversal-invariant gap¬ 
ping of the interface. The bosonic case allows two distinct 
possibilities: (i) K p = m<j x , q p = (0,2), Xp — (g 0) or 
(ii) K p = mcr z , q p = (2,2), Xp = (g 0) with rn an even 
integer and r = 0 or 1. Either possibility yields a minimal 
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charge excitation of e* = 2/m. By the criterion of Ref. [22] 
the interface may be symmetrically gapped when ^X T K~ X Q 
is even. It follows that the smallest possible value of m is four, 
and that the value of r does not affect the interface’s gappa- 
bility. Hereafter we set r = 0 for concreteness and focus on 
K p = 4a x , which relates to the ‘double T-Pfaffian’ discussed 
above (see Appendix). 

To specify the gap-opening interactions we introduce low- 
energy fields describing a given interface y. Right/left-moving 
QSH electron operators are ip r/l, v = e % ^ R i L >v. We use sub¬ 
scripts + and — to denote fields from the adjacent upper and 
lower plates. Operators a± rV = e^ a± ^ and d± iV s e l( ^ d± ^ 
then respectively create self-bosonic charge-e /2 and neutral 
excitations with time-reversal properties a±^ y a±^ y and 
d±, y d± y . These quasiparticles exhibit mutual statistics 
e Z7r / 2 , implying that a± y and d±, y represent local bosons. In¬ 
teractions 


WwM 2 (a_a + ) 4 + H.c. ~ cos AO c ( 2 ) 

(d [_+ H.c. ~ cos 4# s (3) 

^ala + ^ H- H.c. ~ cos 40 nj (4) 

are therefore physical. (We suppress //-dependence whenever 
unneeded.) The fields # c?s?n defined above obey trivial com¬ 
mutation relations and can therefore be simultaneously pinned 
to gap the interfaces. Moreover, the interactions preserve 
both T and charge conservation. Thus uniformly condensing 
(e l o c , s ,n ) ^4 q res p ec t s a n symmetries. 

Determining the resulting surface topological order re¬ 
quires identifying the deconfined anyons. The plates carry 
16 quasiparticles built from combinations of a and d. What 
happens when we drag a quasiparticle from one plate to an¬ 
other as in Fig. 3(b)? Consider first dragging an a charge¬ 
e/2 anyon. Since fractional excitations cannot directly cross 
between plates, this process leaves a dipole described by 
a y- a l + ~ e%dn at ^e interface as the figure illustrates. How¬ 
ever, the condensate {e l6n ) readily absorbs the dipole—which 
is effectively invisible—negating any energy cost. The a 
quasiparticle thus propagates freely across the surface. 

Similarly dragging a neutral d anyon between plates does 
not simply leave behind a d y _ dt dipole since this object is 
uncondensed. To specify its fate we define a neutral fermion 

*pR/L = ^R/LO 2 • (5) 

The condensates identify PpR and Pp\\ we therefore refer to 
both as simply When d crosses an interface it leaves the 
combination d y _d^ y+ ip —which is condensed—and turns into 
a different anyon corresponding to d augmented by the neutral 
fermion pj. Thus quasiparticles d given by 


d = 



even plates 
odd plates 


( 6 ) 


TABLE I. Topological data for the fundamental anyons d and a in 
the symmetrically gapped weak topological insulator surface. 
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also propagate freely across the surface. Remarkably, transla¬ 
tions act nontrivially on these anyons: 

TydT - 1 = #. (7) 

This property, which manifests an even-odd effect, crucially 
distinguishes the symmetrically gapped WTI surface and the 
topological order formed by individual plates. Table I sum¬ 
marizes the topological data for the surface. 

Breaking either translation or time-reversal invariance al¬ 
lows the topological order to appear in strict 2D. Absent the 
former, the gapped WTI surface becomes equivalent to a sin¬ 
gle plate, i.e., the topological order is faithfully described by 
K p = 4a x , q p = (0,2 ),Xp = (0,0). If instead we weakly 
break time reversal the surface should be described by a gen¬ 
erally different K-matrix and charge vector (that again yield 
cr xy = 0 due to the energy gap), as well as a translation matrix 
M y . The latter must obey M y = 1, preserve the K-matrix, 
and encode Eq. (7) while acting trivially on a quasiparticles 
and electrons. Trivial action means that a translated excitation 
at most acquires a local boson that transforms trivially under 
all symmetries. For example, the operation 

MyipnMy 1 = aVjj ~ ip R (v4'0La 4 V’]?V’i) (8) 

multiplies f>R by the local boson in parenthesis. Since this 
boson is odd under time reversal, Eq. ( 8 ) legitimately imple¬ 
ments translations only in the absence of time-reversal sym¬ 
metry. One can indeed view the T-broken surface as a 2D 
system described by [23] 
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with translations implemented in precisely this way. Notice 
that if we ignore M y Eq. (9) could equally well describe a 
time-reversal-invariant 2D phase. However, enforcing transla¬ 
tion symmetry through M y violates time reversal; see Eq. ( 8 ). 
In fact generalizing such 2D realizations to implement time- 
reversal and translation symmetry as in the WTI surface is 
impossible following our earlier general arguments (Fig. 2). 

It is illuminating to discuss the gapped WTI when the sys¬ 
tem consists of a finite stack of N QSH layers. The surface 
is then quasi-ID and hence technically cannot sustain the re¬ 
quired topological order. Indeed, this case reveals a subtlety 
regarding time reversal and the possibility of weak symmetry 
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breaking [24]. As a primer consider Fig. 3(c) where a cylindri¬ 
cal plate ‘wraps around’ a single QSH edge, leaving a gapless 
helical region of length Ax. The QSH/plate interface is iden¬ 
tical to that considered above, and the same interactions (2)- 
(4) can open a gap—ostensibly without breaking symmetries. 
Furthermore the circular edge of the plate can ostensibly also 
be symmetrically gapped (either by finite-size effects or inter¬ 
actions [22]). Interestingly, symmetry-breaking must never¬ 
theless occur [24]: A right-moving electron from the gapless 
QSH edge cannot penetrate into the adjacent gapped segments 
and must reflect into an opposite-spin left-mover, indicating 
spontaneous time-reversal symmetry breaking. 

Using Eq. (3) one can express the magnetization at the 
gapless region’s endpoints as (^ l ^r) ~ e l26s ((d}_d+) 2 ). 
Three cases exist: (i) When interactions gap the plate’s cir¬ 
cular edges the expectation value on the right side must be 
circumference-independent, (ii) When they are gapped by 
finite-size effects, (d}_d+) and thus the magnetization decay 
as a power-law in the cylinder circumference L. (in) Finally, 
when the entire QSH edge is gapped (Ax —)► 0) the circular 
edges are simply absent and ( d)_d+) ~ e~ L ^ with a length £ 
set by the plate’s bulk quasiparticle gap. This last case corre¬ 
sponds to the setup examined in Ref. [24] . Consider next the 
N = 2 generalization of Fig. 3(c) where plates arranged into a 
cylinder gap two QSH layers. The above argument for sponta¬ 
neous time-reversal symmetry breaking no longer holds since 
a right-moving QSH electron from one layer can backscatter 
into the other—which T symmetry readily permits. 

These two examples signify an even-odd effect. For N lay¬ 
ers with periodic boundary conditions along y the local mag¬ 
netization at an interface analogously reads 

~ (^( d -,y d +,v) ^ ' ( 10 ) 

A finite expectation value generically arises if a d quasiparti¬ 
cle from just above the interface can propagate intact to the 
bottom side of the interface. (Direct tunneling is disallowed, 
so the quasiparticle must take the Tong way’ around.) The is¬ 
sue is subtle since d acquires a neutral fermion ^ when cross¬ 
ing an interface; recall Fig. 3(b). With even N the initial d 
ends up dressed by ^ when it reaches the bottom of the inter¬ 
face, and the magnetization thereby vanishes. By contrast, for 
odd N the d quasiparticle boldly arrives undressed yielding 
a finite value. If the entire surface is gapped this expectation 
value decays exponentially with N, while with adjacent gap¬ 
less modes [similar to Fig. 3(c)] a power-law emerges. 

The fact that time-reversal symmetry must be (weakly) bro¬ 
ken in an odd-layer system follows from very general consid¬ 
erations and thus comprises a useful consistency check. The 
surface-state spectrum for an odd number of non-interacting 
QSH layers is gapless with an energy splitting A E oc 1/N to 
the next band. Thus for energies below A E the system maps 
to a single 2D QSH state. Now consider such a finite WTI and 
assume that interactions gap the surface only for x > 0. Elec¬ 
trons from the lowest band must backscatter with unit proba¬ 
bility upon approaching the gapped region; since A E oc 1/N 


the magnetization must also scale as a power-law (or slower) 
to ensure this outcome. This is exactly what we found. 

The non-trivial action of translation symmetry on d anyons 
yields interesting consequences for lattice defects. In a WTI 
screw dislocations terminating at position x 0 on the surface 
(as in Fig. 4) bind a helical QSH edge state that penetrates 
into the bulk [25]. When interactions gap the WTI bound¬ 
ary, electrons from the bulk helical modes must backscat¬ 
ter at the surface. Such a defect thus locally violates time- 
reversal symmetry—yet another manifestation of weak sym¬ 
metry breaking. The impact on surface anyons is even more 
striking: When d encircles the termination point as sketched 
in Fig. 4 it changes any on type and acquires a neutral fermion. 
This suggests that the dislocation forms an extrinsic non- 
Abelian defect that traps a nontrivial zero mode (similar ef¬ 
fects arise in [26-30]). To confirm this, one may formally as¬ 
sociate the QSH edge fields 'ipR/L that enter the bulk with the 
surface modes at x > xo and add a two-particle backscattering 
term (^ r ^l) 2j tH.c. to capture the spontaneous time-reversal 
symmetry breaking. The Appendix derives the following ef¬ 
fective Hamiltonian density that describes the defect, 

U = AQ(x 0 - x)xjj R ^L + uQ(x - xo){^ R tjj L ) 2 + H.C., 

(ID 

with ipR/L defined in Eq. (5). (Note however that at x > xq 
we no longer have ^r ~ ?/>£.) The A and u terms respectively 
arise from Eq. (2) and the two-particle backscattering upon 
taking into account condensates involving the plates. Refer¬ 
ences [31, 32] analyzed precisely Eq. (11) and showed that 
the defect hosts a Z 4 parafermion zero mode. The ‘Z 4 -ness’ 
reflects the two possible values for the spontaneous magneti¬ 
zation and (neutral) fermion parity. 



FIG. 4. (color online) Weak topological insulator surface with a 
screw dislocation. Upon encircling the dislocation d anyons acquire 
a neutral fermion, indicating a zero mode bound to the defect. 

To conclude, we explored strongly interacting WTI surfaces 
using a quasi-ID formulation that permits full analytical con¬ 
trol. We found that the surface can become gapped by en¬ 
tering an Abelian topologically ordered state with several un¬ 
usual features. First, symmetries act on quasiparticles in a 
manner forbidden in purely 2D setups. Second, an interesting 
even-odd effect previously known for non-interacting elec¬ 
trons persists in the topologically ordered surface: For a WTI 
composed of an odd number of QSH systems, ‘weak sym¬ 
metry breaking’ leads to a magnetization exponentially small 
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in the number of layers. Third, lattice defects in the Abelian 
topologically ordered surface exhibit a non-Abelian structure, 
which may be viewed as a manifestation of the anomalous 
symmetry properties of the quasiparticles. We expect such 
features to persist quite generally in weak topological phases 
assembled from 2D symmetry-protected topological states. 
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Derivation of Defect Hamiltonian 

We now derive the Hamiltonian in Eq. (11) that describes 
the dislocation depicted in Fig. 4. Continuing to denote the 
defect position by xo, the gap-opening interactions for x < xo 
are given by 

6c = (Ml) 2 (a-a+) 4 + H.c. 

6 S = (v>]^l) ( [ d -d+) + H.c. ( x < x 0 ) (12) 

d a = (aia+) + H.c., 

which are precisely the terms invoked in the main text to gen¬ 
erate symmetric topological order. At x > xq the upper and 
lower plates fuse together without the aid of an intervening 
QSH edge mode—which detours into the bulk at xo as the 
figure illustrates. As in the main text we use x > xo to denote 
both the QSH edge modes that enter the bulk and plate fields 
at the surface. The following perturbations describe the fusion 
of the plates in this region, 

o d =(did+y+ h.c. 

O a = (a f _a+) +H.c., (x>x 0 ). (13) 

When relevant the above terms catalyze condensation of 
(d}_d+) and {a}_a+), allowing the anyons to seamlessly pass 
between plates as desired. It remains to specify the fate of the 
QSH edge fields ^ r/l that bleed into the bulk. We will add a 
perturbation (assumed relevant) 

Of = (^IjV’l) + H.C., (X > :ro) (14) 
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for those modes; this catalyzes a magnetization 7^ 0 

without changing the symmetries of the problem since the 
defect spontaneously breaks time-reversal symmetry as dis¬ 
cussed earlier. 

Notice that O a is present on both sides of the defect. We 
can therefore set a = a_ ^ a + in all remaining terms. Fur¬ 
thermore, O s at x < xo is simply a product of terms present 
in Od and O^ for x > x$. Minimizing the energy therefore 
requires fixing to a uniform constant every¬ 

where. (Otherwise there will be an energy cost when the ex¬ 
pectation value ‘twists’ at xo-) It follows that d}_d+ is slaved 
to the magnetization in the region x > xq. 

All the interesting defect physics has now been distilled into 
the perturbations O c and O^. In terms of neutral fermions 
i’R/L = i>R/La 2 we have 

(Ml) 2 + H.c. (15) 

and 

O* ~ (i>y L y + H.C. (16) 

Focusing on these crucial pieces, the defect Hamiltonian den¬ 
sity becomes 

H' = A0(x o - zXVw/’z,) 2 + uQ(x - x 0 XV’JjV’i,) 2 + H c - 

The first and last terms respectively favor condensates with 
$r^l) 0 and 7^ 0- The sign of the latter is ar¬ 

bitrary but, importantly, the former is not: The paired neutral- 
fermion condensate Josephson couples to the remaining 2D 
surface and hence its phase is chosen spontaneously but glob¬ 
ally. For the ID subsystem describing the zero modes it thus 
acts like an external superconductor with a fixed phase, and 
one should replace ^ We then arrive pre¬ 

cisely at the Hamiltonian density quoted in Eq. (11). 


Relationship Between Symmetrically Gapped Strong and Weak 
Topological Insulator Surfaces 

Here we comment further on the connections between 
the symmetric Abelian topological orders for gapped WTI 
surfaces and the known non-Abelian topological orders for 
gapped STI surfaces. For the latter, two phases unrelated 
by simple condensation transitions are possible [6-9] : the T- 
Pfaffian noted in the main text and the somewhat more elab¬ 
orate Pfaffian-antisemion phase. Four possibilities for gap¬ 
ping the WTI surface thus immediately arise—each of the 
two Dirac cone can be gapped independently into either of 
these phases. As we argue below, we expect that both of the 
minimal Abelian topological orders that can arise on the WTI 


surface—corresponding to plates with K = 4a x and 4a z — 
may be accessed from these non-Abelian states via condensa¬ 
tion transitions. Which Abelian states arise depends on which 
of the four possible parent non-Abelian phases one considers. 
In particular, for the simplest non-Abelian ‘double T-Pfaffian’ 
state, obtaining topological order corresponding to 4a z can be 
ruled out leaving 4a x as the only natural possibility. 

To relate the possible STI and WTI states, suppose first that 
the WTI surface can be fully gapped by purely local interac¬ 
tions between 2D topologically ordered plates and QSH edges 
(i.e., acting solely within a given interface). This is precisely 
the situation we encountered (by construction) in the main 
text. In this case the WTI surface inherits only the quasipar¬ 
ticles native to the 2D plates, with no additional topologically 
distinct excitations generated. Continuing to assume the mini¬ 
mal two-dimensional K-matrix for the plates, the quasiparticle 
content of the symmetric topologically ordered WTI surface 
thus corresponds to either K = 4 a x or 4a z as discussed be¬ 
low Eq. (1). 

Suppose next that the symmetrically gapped WTI surface 
instead descended from a topologically ordered parent state— 
e.g., the ‘double T-Pfaffian’—via a simple condensation tran¬ 
sition. From this viewpoint the quasiparticles and their prop¬ 
erties are inherited from the parent phase, and simply corre¬ 
spond to the anyons that remain deconfined after the transi¬ 
tion. Again, no additional anyons appear in the process. 

We can thus rule out potential parent states by identifying 
quasiparticles that appear in the descendant state, but not in 
the putative parent. Symmetric WTI surface topological or¬ 
ders corresponding to 4a x and 4a z both admit quasiparticles 
with charge e/2. In the 4a x case these excitations carry topo¬ 
logical spins e m 5 for integer n, while for 4 a z their topologi¬ 
cal spins are instead given by e ^ 2n+1 ^ f. Regarding the parent 
states, both the T-Pfaffian and Pfaffian-antisemion phases host 
Abelian e/2 excitations with topological spin =t i; the Pfaffian- 
antisemion (which admits more anyon types) hosts additional 
e/2 quasiparticles with spins ±1. Both phases also harbor 
non-Abelian e/4 excitations, which carry spin 1 in the T- 
Pfaffian but spin e 27r / 4 in the Pfaffian-antisemion [8]. Thus in 
either the ‘double T-Pfaffian’ or ‘double Pfaffian-antisemion’, 
all quasiparticles with half-integer charge have topological 
spins of the form e m ^, consistent with 4a x but not with 4a z . 
In contrast, a phase where the two Dirac cones are gapped 
differently—i.e., one in the T-Pfaffian and one in the Pfaffian- 
antisemion—exhibits e/2 excitations with topological spins 
e m7r / 4 . Hence this phase could potentially give way to 4 a z 
upon condensing appropriate combinations of quasiparticles. 

We note that the analysis presented in main part of the paper 
may be interpreted as explicitly carrying out a condensation 
starting form a symmetric parent state that is closely related 
to the ‘double T-Pfaffian’. In this modified parent state, the 
two copies of T-Pfaffian topological order are not indepen¬ 
dently symmetric under time reversal and translations, but in¬ 
stead related by these symmetries; see Fig. 1(b) and Ref. [11]. 
This approach has the technical advantage that it permits a 
controlled analytic treatment in terms of quasi-ID Hamiltoni¬ 
ans and bosonization. Apart from this, we do not expect any 
significant difference between either method. 


